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' Abstract. We derive, for a bistochastic strictly contractive quantum channel on 

a matrix algebra, a relation between the contraction rate and the rate of entropy 
production. We also sketch some applications of our result to the statistical physics of 
C*") ■ irreversible processes and to quantum information processing. 
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Let A be the algebra of observables (say, a C*-algebra with identity), associated 
with a quantum-mechanical system E. A general evolution of S is described, in the 
Heisenberg picture, by a map T : A — » A which is (i) unital: T(l) = 1, and (ii) 
completely positive: for any nonnegative integer n, the map T®id : A® M. n — > A®M. n , 
where Ai n is the algebra of n x n complex matrices, sends positive operators to positive 
>• . operators Q. We shall henceforth refer to such maps as (quantum) channels 0. If 
^ . S is an iV-level system, then its algebra of observables is isomorphic to AAn, which 
is exclusively the case we shall consider in this letter. A celebrated result of Kraus 
[ lj then says that, for any channel T, there exists a collecton of at most N 2 operators 
Vi G M.Ni which we shall call the Kraus operators associated with T, such that (i) 
T ( A ) = Ei VfAVi, and (ii) £\ V*Vi = 1. It is now easy to see that T(A*) = T(A)* for 
all A G M.Ni i-e., any channel maps Hermitian operators to Hermitian operators. 

Given a channel T, the corresponding Schrodinger-picture channel T is defined via 
the duality 

tr[f(A)B] = tv[AT(B)}, 

whence it follows that T is a completely positive map which preserves the trace, i.e., 
trT(A) = tr A for all A E Mn- in other words, T maps the set of N x N 
density matrices into itself. Furthermore, in terms of the Kraus operators Vi, we have 
T{A) = E 4 ViAV*, so that f(A*)= f(A)* as well. 

The set {T n } ng N is a discrete-time quantum dynamical semigroup generated by T, 
i.e., T n T m = T n+m and we take T° = id (the identity channel). It is easy to show that, 
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for any channel T, there exists at least one density operator p such that T(p) = p [H|. A 
question of clear physical importance is to determine whether the dynamics generated by 
T is relaxing i.e., whether there exists a density operator p such that, for any density 
operator a, the orbit {T"(er)} converges to p in the trace norm ||A||i:= tr (A* A) 1 / 2 . 

One way to show that a dynamics is relaxing relies on the so-called Liapunov's 
direct method ||. Let X be a. compact separable space, and let <p : X — > ^ be a 
continuous map, such that 

(i) ip has a unique fixed point xq G Af, and 

(ii) there exists a Liapunov function for 99, i.e., a continuous functional 5 on ^ such 
that, for all x E X , S[ip(x)] > where equality holds if and only if x = x . 

Then, for any x G X, the sequence {<p n (x)} converges to x . 

Let T be a bistochastic channel, i.e., one for which T(l) = T(lL) = 1. If we treat 
.Mat as a Hilbert space with the Hilbert-Schmidt inner product, (A, B):=tr (A*B), then 
an easy calculation shows that the Schrodinger-picture channel T is precisely the adjoint 
of T with respect to (•,•), i.e., (A, T(B)) = (f(A),B) for all A,B e M n . The composite 
map T o T (which we shall write henceforth as TT) is also a bistochastic channel, which 
is, furthermore, a Hermitian operator with respect to (•,•)■ In ||, Streater proved the 
following result. 

Theorem 1 Let T : M.n — ► M-n be a bistochastic channel. Suppose that T is ergodic 
with a spectral gap 7 G [0, 1), i.e., (i) up to a scalar multiple, the identity matrix 11 is 
the only fixed point of T in all of M.n, and (ii) the spectrum of TT is contained in the 
set [0, 1 — 7] U {1}. Then, for any a G TV, we have 

S[f{a)\-S{a)>l\\a-N-Hf 2 , (1) 

where S{p):= — tr (oTncr) is the von Neumann entropy of a and ||v4||2: = [tr (v4*v4)] 1//2 is 
the Hilbert-Schmidt norm of A. 

In other words, if a bistochastic channel T is ergodic, then the dynamics generated 
by T is relaxing by Liapunov's theoremf]]. Furthermore, the relaxation process is 
accompanied by entropy production at a rate controlled by the spectral gap. 

Now we have an interesting "inverse" problem. Consider a bistochastic channel 
T on M.n with T strictly contractive 0. That is, T is uniformly continuous on 
T>n (in the trace- norm topology) with Lipschitz constant C G [0, 1): for any pair 
a, a' G V N , we have ||T(er) — T(cr') ||i < C\\a — er'||i. Then by the contraction mapping 
principle J7J, iV" 1 ! is the only density matrix left invariant by T, and furthermore 
||T(cr) — 1 H || 1 — ► as n — ► 00 for any cr e V N , i.e., the dynamics generated by T is 
relaxing. The question is, does the entropy-gain estimate (|1|) hold, and, if so, how does 
the spectral gap 7 depend on the contraction rate C? 

f Endowing the set 2? at with the trace- norm topology takes care of all the continuity requirements 
imposed by Liapunov's theorem. 
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This problem was motivated in the first place by the following observation. In the 
case of M.2-, the action of a bistochastic strictly contractive channel T can be given a 
direct geometric interpretation. Recall that the density matrices in are in a one- 
to-one correspondence with the points of the closed unit ball in R 3 . Then the image 
of T>2 under a strictly contractive channel T with contraction rate C will be contained 
inside the closed ball of radius C, centered at the origin 0, i.e., the image of V 2 under 
T will consist only of mixed states. This geometric illustration suggests that the rate of 
entropy increase under T must be related to the contraction rate. Now even though in 
the case of M. n with N > 3 we no longer have such a convenient geometric illustration, 
nevertheless it seems plausible that the rate of entropy production under a bistochastic 
strictly contractive channel would still be controlled by the contraction rate. 

Indeed it turns out that the contraction rate is related to the rate of entropy 
production, as stated in the following theorem. 

Theorem 2 Let T be a bistochastic channel on M.N-, such that T is strictly contractive 
with rate C . Then T is ergodic with spectral gap 7 > 1 — C, so that, for any a G Vn, 
we have 

Slfia^-Si^y^Wa-N-HWl. (2) 

Proof. We first prove that T is ergodic. As we noted before, T and T are adjoints 
of each other with respect to the Hilbert-Schmidt inner product. Using the Kadison- 
Schwarz inequality M 

<$>(A*A) > $(A)*Q(A) 
for any channel $ on a C*-algebra A, as well as the fact that 

trT(A) = tr[f(lL)A] = tr A 

for any i 6 M;v, we find that 

\\T(A)\\l = tr [T{A)*T{A)\ < tr[T(A*A)} = tr (A* A) = \\A\g 

and the same goes for T. That is, both T and T are contractions on A^jy (in the 
Hilbert-Schmidt norm), hence their fixed-point sets coincide ||. 

By hypothesis, T leaves invariant the density matrix iV -1 !!, which is invertible. In 



this case a theorem of Fannes, Nachtergaele, and Werner |TI| |Tl]] says that T(X) = X 
if and only if VJX = XVi for all V$, where Vi are the Kraus operators associated with 
T. It was shown in |J that if T is strictly contractive, then the set of all X such that 
ViX = XVi for all Vj consists precisely of multiples of the identity matrix. We see, 
therefore, that T(X) = X if and only if X = %1 for some \ G C, whence it follows that 
T{X) = X if and only if X is a multiple of 1. This proves ergodicity of T. 

Our next task is to establish the spectral gap estimate 7 > 1 — C. Let X be a 
Hermitian operator with trX = 0. In that case we can find a density operator p and 
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a sufficiently small number e > such that a:=p + eX is still a density operator ||||. 
Then 

llTWHi = (lA)ll^) - T(p)h < (C/e)\\a - p\\ x = CftX^. (3) 

Because one is a simple eigenvalue of both T and T, it is also a simple eigenvalue of TT. 
Hence 1— 7 (which we may as well assume to belong to the spectrum of TT) is the largest 
eigenvalue of the restriction of TT to traceless matrices. Let Y be the corresponding 
eigenvector. Without loss of generality we may choose Y to be Hermitian||. Then, using 
Eq. (HD and the fact that ||$(A)||i < ||^4||i for any trace-preserving completely positive 
map $ jn|], we may write 

(i-7)imii = l|TT(y)|| 1 <||f(y)|| 1 <c'||y|| 1 , 

which yields the desired spectral gap estimate. The entropy gain bound @ now follows 
from Theorem [l]. I 



Remark. Eq. (|3]) can also be proved using the following finite-dimensional 

specialization of a general result due to Ruskai |13[]. If T : M. N — > M. N is a channel, 
then 

\\f{A)\\i 1 

A=A*;trA=0 ||^4||l ^ i/j,<f>eC N ;{il}\(f>)=0 



sup }L JmT ± = 2 SUP IWXV'l-I^X^DIIi- (4) 



Because T is strictly contractive, the right-hand side of (f|) is bounded from above by 
C, and d||) follows. The supremum on the left-hand side of (§) is the "Dobrushin con- 
traction coefficient," studied extensively by Lesniewski and Ruskai 0] in connection 
with the contraction of monotone Riemannian metrics on quantum state spaces under 
(duals of) quantum channels. □ 



Note that in some cases the sharper estimate 

S[f(a)]-S(a)>^-^\\a-N-H\\l (5) 

may be shown to hold. Consider, for instance, the case T = T, so that the eigenvalues 
of T are all real. Let Ai, . . . , Xl, L = N 2 — 1, be the eigenvalues of T that are distinct 
from unity. Then we claim that max, \Xj\ < C, which can be proved via reductio ad 
absurdum. Suppose that there exists some X (which we may take to be Hermitian) 
with txX = such that T(X) = XX with |A| > C. We may use the same trick 
as in the proof above to show that there exist two density operators, a and p, such 
that ||T(cr) — r(p)||i > C||er — p||i, which would contradict the strict contractivity of T. 

\ This may be seen as a simple consequence of the following fact . The set of all invertible 
N x N density matrices is a smooth manifold, where the tangent space at any p 6 T>™ can be naturally 
identified with the set of N x N traceless Hermitian matrices. 

§ Recall that 1 - 7 is real, and TT(A*) = [TT(A)]* for all A, which implies that Y + Y* is also an 
eigenvector of TT with the same eigenvalue. 
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Because TT = T 2 , we have 1— 7 = (max.,- |Aj|) 2 < C 2 , which confirms (|5]). Furthermore, 
using a theorem of King and Ruskai the bound ([5]) can be established for all 
bistochastic strictly contractive channels on M.2-, as well as for tensor products of such 
channels. 

The proof of this last assertion goes as follows. Let T be a channel on .M2 such that 
T is strictly contractive. Then T is ergodic, the proof of which can be taken verbatim 
from the proof of Theorem ||. It is left to show that 1 — 7 < C 2 . Any density operator 
in M.2 can be written as 

where the o~j are the Pauli matrices 

/ -i \ ( 1 

-1 = 1 , n 1, ^=^. j, - 3 = ^ _ x 

and the real numbers rj satisfy the condition r\ + r\ + r\ < 1 (this is precisely the one- 
to-one correspondence between T>2 and the closed unit ball in IR 3 ). The King- Ruskai 
theorem |T5j asserts that, for any bistochastic channel T on M.2, there exist unitaries 




U, V and real numbers £j, 1 < j ' < 3, with |^-| < 1 such that, for any p e V 2 , 

f(p) = U[f dias (VpV*)]U*, 
where the action of the map Tdi ag on the density operator (g) is given by 



T, 



diag 



It can be easily shown that if T is strictly contractive, then C = max^ |£j 



The parameters £j are determined as follows [15j. Consider the orthonormal basis 
of M.2, generated by 11 and the Pauli matrices, with respect to which T and T can be 
written as 4 x 4 matrices in the block-diagonal form 




T= I . I , T 



1 

M* 

where M* denotes the transpose of M0J Then the absolute values of the parameters 
£j are precisely the singular values of M. This implies that C = \\M\\, where ||M|| 
denotes the operator norm (the largest singular value) of M. Consequently we get 
1 — 7 = ||M*M|| = ||M|| 2 = C 2 , which yields the entropy gain estimate (|5]). The 
argument for tensor products of bistochastic channels 011M2 runs along similar lines. 

If The 3x3 matrix M must be real because both T and T map Hcrmitian operators to Hermitian 
operators. 
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We remark that the results reported in this letter are consistent with the following 
theorem ||. Let T be a channel on A4n with the property that T has a unique fixed 
point p G T>n in all of A4n- Let Xj be the eigenvalues of T distinct from one, and let 
/t:=maxj |Aj|. Then there exist a polynomial p and an iV-dependent constant K such 
that, for any a G £V, 

||T>)-p||i<ii>(nK. (7) 

This shows that the dynamics generated by T is relaxing, ||T n (cr) — p||i — > as n —>■ oo, 
and that the rate of convergence is controlled essentially by the eigenvalue of T with the 
second largest modulus. Now, if T is a bistochastic channel with T strictly contractive, 
then it follows from Theorem |2| that k < C 1 ^ 2 . But by virtue of strict contractivity we 
have 

2C n / 2 fAT- n 

||T"(a) — iV = ||T n (a) — T n (iV l)||i < C n \\a - N' 1 ^ < — 



N 

which has the form of (|7|). To obtain the last inequality we used the fact C < C 1 ^ 2 
for < C < 1, as well as the fact that the set T>n is compact and convex, so that the 
convex functional ||er — A^ -1 !^ attains its supremum on an extreme point of T>n, i.e., 
on a pure state. In turn, a routine calculation shows |6j that, for any pure state a, 
\\a - N- 1 ^ = 2(JV- 1)/N. 

Finally, a few comments are in order as to how our results come to bear upon (a) the 
statistical physics of irreversible processes 0] and (b) quantum information processing 

In the first setting we are interested in a concise mathematical description of 
the approach to equilibrium. We consider a quantum-mechanical system with the 
Hilbert space TC, dim H < oo. We will use A to denote the algebra B(TC) of all 
bounded linear operators on TC. Let H be the (possibly unbounded) Hamiltonian 
of the system, with the additional requirement that the operator e~^ H is trace-class 
for all positive real (3. This means that the Gibbs state exists for all positive inverse 
temperatures, and that each eigenvalue of H has finite multiplicity. Let us write down 
the spectral decomposition H = J2e Pe> where E are the eigenvalues of H and Pe are 
the corresponding eigenprojections. Then he'-= dim PeTC < oo for all E. Suppose we 
are given a channel T : A — > A such that, for each E, the restriction Te of T to the 
algebra Ae'—PeAPe is bistochastic, i.e., T(Pe) = T(Pe) = Pe, so that T(Ae) Q Ae, 
and strictly contractive. Then Theorem |^ says that, for each E, Eq. (Q) holds with some 
contraction rate Ce- The Hilbert spaces PeTC can be thought of as the energy levels 
of H, and the states PeI^e as the corresponding microcanonical states. If in addition 
we can show that C:= sup B Ce < 1, then any normal state of the entire system will 
converge toward the Gibbs state with the rate C. 

In the framework of quantum information processing we are concerned mostly with 
systems that have finite-dimensional Hilbert spaces. We may envision either a quantum 
register that is used for storing information, or a quantum computer that processes 
information. In both of these situations we are interested in stabilizing the information 
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against the effects of decoherence. The decoherence mechanism is modelled by a channel. 
It was shown in || that the set of strictly contractive channels on an algebra A is dense in 
the set of all channels on A (here we are talking about the Schrodinger-picture channels). 
The same result also holds for bistochastic strictly contractive channels, which are dense 
in the set of all bistochastic channels |J. For such channels Theorem [2] may be used to 
estimate the rate at which entropy is produced in the register or in the computer. Such 
estimates generally serve as a measure of efficiency of error- correct ion procedures [fL7f 
(note that it was shown in || that errors modelled by strictly contractive channels can 
be corrected only approximately). 
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